Monopole oscillations in light nuclei with a molecular dynamics approach 
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Collective monopole vibrations are studied in the framework of an antisymmetrized version of 
molecular dynamics as a function of the vibration amplitude. The giant monopole resonance energy 
in 40 Ca is sensitive to the incompressibility of the effective interaction, in good agreement with 
complete time-dependent Hartree-Fock calculations. The collective response of 12 C, 16 0, and 24 Mg 
is also studied. For these lighter nuclei that have an important contribution of an a-clustered 
component, different frequencies are observed, corresponding to two different types of vibrations 
associated with breathing and moving of the underlying clusters. Possible connections with direct 
breakup into alpha clusters at high excitation energy are discussed. 

PACS numbers: 24.30.Cz,21.65.-f,24.30.Gd 



I. INTRODUCTION 

Collective vibrations of nuclei are interesting phenom- 
ena in quantum many-body dynamics. The constituent 
nucleons coherently participate in these vibrations and 
the motion can be intuitively understood in classical 
terms as a global vibration of an excited liquid drop. In 
a quantum mechanical sense, the dynamics is described 
by a coherent superposition of many particle-hole exci- 
tations on top of the ground-state configuration. These 
collective vibrations have been experimentally observed 
for thirty years as broad resonances (giant resonances) in 
the high-energy spectrum of nuclei pH4[. The existence 
of such resonances has been reported for many different 
medium-heavy and heavy nuclei in the isotopic table, and 
their characteristics have been successfully reproduced in 
mean-field-based approaches Q. 

These theoretical studies have shown that collective 
modes are privileged probes of the effective interaction, 
and as such they can bear important information on nu- 
clear matter properties. In particular, the isoscalar gi- 
ant monopole resonance or the breathing mode of nuclei 
is connected to the incompressibility of nuclear matter, 
even if this connection is far from being direct @, Q • 

Experimental data on the isoscalar monopole response 
of nuclei as light as 12 C or 16 (8rJIT| show a fragmented 
strength distribution, which is qualitatively understood 
as a global weakening of the collectivity but is not quan- 
titatively understood. It is well known that correlations 
play an important role in the structure of these light nu- 
clei and the application of standard Hartree-Fock (HF) or 
random phase approximation (RPA) calculations to their 
collective dynamics may be questionable. This is partic- 
ularly true in nuclei presenting clustered structures, a 
phenomenon that is expected both in excited states and 
in the ground state of numerous light and medium-heavy 
N = Z nuclei [HI, [l3|. In such situations, one can ex- 
pect that collective vibrations may imply collective mo- 
tions of the underlying clusters. Indeed, it has already 
been reported in the literature that collective strengths 
to low-lying states are affected by the cluster degrees of 



freedom [3 (3. 

The giant monopole resonance (GMR) corresponds to 
relatively small amplitude of a collective monopole exci- 
tation. Another interesting regime of nuclear dynamics is 
given by large-amplitude collective motion with a signif- 
icant compression of the system followed by expansion. 
At excitations of the order of 3 AMeV or higher, it leads 
to the observable phenomenon of nuclear fragmentation 

Here the situation is qualitatively similar to the excita- 
tion of giant resonances: When a heavy nucleus is com- 
pressed significantly, the dynamics of compression and 
expansion appears to be governed by global collective 
properties, namely, the nuclear matter equation of state, 
and appears to be largely independent of the structure 
properties of the excited nucleus 0, HU . Conversely, if 
the ground state and/or excited states have pronounced 
clustered structure, one may expect an interplay between 
these structures and the collective dynamics and frag- 
mentation pattern. A systematic study of the fragmen- 
tation of light clustered nuclei has never been done, but 
some indications exist as a preferential a-particle emis- 
sion from a-structured nuclei [22l - |25| . 

In this respect, models able to describe both the small- 
and large-amplitude collective regime are useful tools 
for understanding collective nuclear dynamics. Interest- 
ing studies have been performed in the context of time- 
dependent Hartree-Fock (TDHF) calculations [26l4a3| 
and its semiclassical counterpart 35 38] and a cluster- 
based time-dependent approach [39144 1| . 

Antisymmetrized versions of molecular dynamics mod- 
els, antisymmetrized molecular dynamics (AMD) [42l.l43l| 
and fermionic molecular dynamics (FMD) [44j |. are pow- 
erful approaches for describing both nuclear structure 
and reactions. They are microscopic models based on 
nucleonic degrees of freedom. The wave function is mod- 
eled as a fully antisymmetrized product of Gaussian wave 
packets. The choice of the nonorthogonal Gaussian wave 
packets is a restriction of general Slater determinants, but 
it has an advantage of allowing spontaneous breaking of 
all spatial symmetries. Because of that, these models are 
particularly suitable to describe cluster and molecular 
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states. After projection over appropriate quantum num- 
bers, they have been successfully applied to describe the 
structure change from shell-modcl-like structure to clus- 
tering structure [44M46| . They also have been applied 
to nuclear reactions from fusion (47j to multifragmenta- 
tion reactions (U EH, EH . For small-amplitude collective 
modes, the isovector dipole oscillation has been studied 
with AMD [49| but almost no application to the isoscalar 
monopole modes has been performed to our knowledge. 

In this paper, we use the simplest version of the model 
employing a single Slater determinant. This approxi- 
mation is not sufficient to realistically describe nuclear 
spectroscopy, but it allows us to understand qualitatively 
the effect of clustering on the collective dynamics. The 
use of a fully time-dependent approach allows us in prin- 
ciple to study consistently the motion from the small- 
amplitude RPA regime to the large-amplitude fragmenta- 
tion regime. In this first application, two-body collisions 
and quantum branching, which are necessary to explain 
experimental data of heavy-ion collisions around Fermi 
energies [H, El, El, H3| , ar e not implemented, but their 
inclusion is straightforward. 

The main difference between AMD and FMD concerns 
the width degree of freedom. In FMD, the centroids and 
width of Gaussians are both complex time-dependent 
variational parameters whereas, in AMD, the width is 
an external real number common to all the Gaussians. 
The isoscalar monopole oscillation is supposed to corre- 
spond to a breathing mode, so that the width degree of 
freedom may have a significant importance. By compar- 
ing the results between AMD and FMD, we would also 
like to clarify what degrees of freedom are important to 
properly describe compression modes. 



II. FORMALISM OF AMD AND FMD 



where bi (i — I,..., A) are complex time-dependent 
variational parameters. In the FMD formulation, the 
wave packet width ai is an additional complex time- 
dependent variational parameter. In contrast, in the 
AMD model, all wave packets share the same real con- 
stant width a,; = l/2v = const. The spin-isospin states 
(Xa; — P tiP -l, n t, n I) give additional variational pa- 
rameters, which are essential to describe charge-exchange 
reactions and spin excitations. In the present application 
to scalar-isoscalar collective modes, only spin and isospin 
saturated systems will be considered. We will then re- 
strict ourselves to a simple variation in the configuration 
space and the spin-isospin states will be treated as con- 
stants of motion. The equations of motion for AMD and 
FMD as resulting from Eq. (JTJ) can be written as 



dqj T _ d{H) 



(4) 



where (H) is the expectation value of the Hamiltonian H 
taken with the wave function VP, Ci a .j T appears because 
of the non-orthogonality of Gaussian wave packets and is 
given by 



if 



■ln(W|W) 



(5) 



and qi a give the variational parameters for nucleoli i in 
the model, namely, qi a =i,...,i = o,i,bi X ,bi y ,bi Z in FMD 
and <7i CT =i 3 = bi X ,bi V ,bi Z in AMD. For further details, 
see Refs. |H,E1- 



III. DESCRIPTION OF GROUND-STATE 
PROPERTIES 



In the AMD and FMD models, the optimal time depen- 
dence of the iV-body wave function is obtained through 
the time-dependent variational principle 



*(t)\ih£ - H\*(t)) 



dt = 



(1) 



(*(*)!*(*)) 

<Wi)> = *|*(ta)> = 

within the constraint that \^(t)) is given by a Slater de 
tcrminant of the single-particle states -0, 

1 



I*) 



= det 0')] 
1 A\ « 



(2) 



where A is number of nucleons. 

If no extra restriction on the single-particle wave func- 
tions 4>(r) is applied, Eq. ((TJ) gives the TDHF equation 
[5l| . However, in the AMD and FMD models, the single- 
particle wave functions rp(r) are restricted to the Gaus- 
sian form: 



ipi (r) = exp 



-i 1 



Xo 



(3) 



The ground state is obtained by minimizing the energy 
expectation value (H) = with respect to 

the variational parameters of the models. 

The molecular dynamics ansatz Eq. ([3]) represents a 
restriction of the full space for the single-particle wave 
function. This means that a full HF calculation in three- 
dimensional r-space is a generalization of this model, 
though molecular dynamics can be superior to HF if the 
latter is solved on a truncated basis or by imposing the 
conservation of chosen spatial symmetries [52j . As a con- 
sequence, AMD and FMD ground-state energies have to 
be seen as upper limits of the HF ground-state value. 

Table U shows some ground-state properties, namely, 
the ground-state energy -E'g.s. and the matter root-mean- 
square radius R m obtained by AMD and FMD for sev- 
eral selected nuclei. Results with SLy4 f53l . Gogny [54| . 
and the old SIII effective interaction |55( are shown. 
SLy4 and SIII are Skyrme interactions, namely, of 5- 
type, with symmetric nuclear matter incompressibilities 
of Koo = 230 MeV and = 355 McV, respectively. 
Gogny is a finite-range interaction with an incomprcss- 
ibility modulus of = 228 McV. Both Gogny and 
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SLy4 are realistic interactions that have been success- 
fully compared to different experimental data on nuclear 
masses and collective modes. SIII also gives reasonable 
values for nuclear masses along the stability valley, but 
it is known to exhibit a too high incompressibility. The 
results of a three-dimensional HF calculation with SLy4 
using the code of Ref. [52j and the experimental data 
from the Atomic Mass Evaluation [56[ are also shown in 
the table for comparison. 





AMD 


FMD 


HF 


Exp 


SLy4 


SIII 


Gogny 


SLy4 


SIII 


SLy4 




4 He 


Eg. s . 


-26.0 


-25.6 


-28.9 


-26.2 


-25.5 


-26.7 


-28.3 


Rm 


1.73 


1.69 


1.63 


1.73 


1.69 


1.96 




12 C 


-Eg. s . 


-72.2 


-70.6 


-74.9 


-77.4 


-75.2 


-90.7 


-92.2 


Rm 


2.57 


2.53 


2.52 


2.57 


2.51 


2.44 




16 Q 


-Eg.s. 


-121.5 


-120.4 


-125.3 


-127.9 


-125.4 


-128.6 


-127.6 


Rm 


2.63 


2.59 


2.58 


2.63 


2.60 


2.67 




40 Ca 


Eg.s. 


-330.6 


-322.1 


-334.2 


-338.9 


-330.3 


-344.6 


-342.1 


Rm 


3.40 


3.42 


3.36 


3.40 


3.41 


3.39 





TABLE I: Ground-state properties [energy E s . a . (MeV) and 
root-mean-square radius R m (fm)] for 4 He, 12 C, 16 and 40 Ca 
calculated by AMD and FMD with three different effective 
interactions (SLy4, SIII, and Gogny). Calculation results of 
the HF method with SLy4 using the code of Ref. [13] and 
experimental data (EXP.) from the Atomic Mass Evaluation 
[5(| are also shown in the table for comparison. 



AMD ^ g . s . 


SLy4 


SIII 


Gogny 


4 He 


0.19 


0.20 


0.21 


12 C 


0.155 


0.16 


0.16 


16 


0.155 


0.16 


0.16 


40 Ca 


0.125 


0.125 


0.13 



TABLE II: List of the inverse width optimal value fg. s . (fm 2 ) 
for the AMD calculations. 

In AMD calculations, the inverse width v is an external 
parameter optimized for each nucleus to ensure a mini- 
mum (H). The optimal values v g . s . are listed in Table [TT1 
The main difference between AMD and FMD in this con- 
text, therefore, is that (i) AMD takes a common width 
for all the Gaussians whereas FMD can take a differ- 
ent width for each Gaussian and (ii) the imaginary part 
of the width is an extra variational parameter for FMD 
whereas it is fixed to zero in AMD. In actual applications, 
point (i) is the most important since it helps in describ- 
ing neutron skins and neutrons halos [57H59| . However, 
the imaginary parts of the variational parameters turn 
out to be very close to zero in FMD ground states, even 
if they can play a role in the dynamical evolution. 

As was already observed in previous more systematic 
studies [ii], Sll , the overall experimental features of the 



ground-state properties are reproduced by both AMD 
and FMD. The difference between the results of HF and 
those of AMD and FMD is not significant for most of 
the cases, even though the state |^) in AMD and FMD 
does not cover the whole possible range of Slater deter- 
minants. This result shows that the Gaussian ansatz has 
a high degree of flexibility, at least as far as the shown 
observables are concerned. 

The only exception is given by the case of 12 C, which is 
significantly underbound by both AMD and FMD. This 
discrepancy is not due to the limited variational space of 
the spatial wave functions, but rather to the lack of the 
spin-orbit interaction, which has been neglected in the 
present work. In fact, if we switch off the spin-orbit in- 
teraction in the HF calculation, we obtain E g . s , = —74.8 
MeV, which is close to the FMD result. This value is 
slightly higher than the FMD result, which may be due to 
the mirror symmetry in xyz axes imposed on the single- 
particle wave functions in the HF calculation [52|. To 
properly introduce the spin-orbit interaction into FMD 
(and AMD), the variation of the spin wave function to- 
gether with the configuration wave function [44| , as well 
as the angular momentum projection, would be neces- 
sary j57H60| . These features are essential for a quanti- 
tative study of nuclear structure and spectroscopy, but 
they are not expected to play an important role in the 
study of the monopole oscillation. 

Because of the restriction of the accessible configura- 
tion space in molecular dynamics models, the following 
relation is expected among the ground-state energies cal- 
culated by these models if the same effective interaction 
is used: 

E s ,s. (HF) < Eg*. (FMD) < E g . s . (AMD). (6) 

However, this inequality is not systematically observed. 
This is due to the different treatment of the spurious zero- 
point kinetic energy of the center-of-mass motion T zero . 
It is well known that the center-of-mass motion gives rise 
to a spurious contribution in HF-based models, owing to 
the breaking of translational invariance. The first-order 
contribution can be subtracted by replacing the kinetic 
part of the Hamiltonian i?kin = 57^; Sj=i ft 2 with 

1 f A ~\ 2 

(7) 

1 ( A A \ 

where m is the mass of the nucleon. The one-body term 
is readily implemented by multiplying a factor (1 — 1/A) 
to the kinetic energy. The two-body term is numerically 
expensive, and it is neglected both in the standard choice 
of HF [52| and in our FMD calculations. Conversely, in 
AMD, the center-of-mass coordinate can be separated 
from the intrinsic degrees of freedom and there is no 
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spurious coupling between them. Therefore, the contri- 
bution of the spurious center-of-mass motion is readily 
calculated as 



1 (*AM D |(Etl^) 2 |*AMD> _ Zh 2 V 



2mA 



amdI^amd 



2m 



(8) 



and can be exactly subtracted from the AMD ground- 
state energy. This approximate treatment of the center- 
of-mass motion in HF and FMD can be of some im- 
portance in the case of the light nuclei considered here. 
When employing in the AMD calculation the same ap- 
proximate treatment of the center-of-mass motion that 
is currently used in FMD and HF, we have verified that 
indeed E g . s .(HF) < £ g . s .(FMD) < £ g . s .(AMD) is always 
found. 



IV. DESCRIPTION OF THE ISOSCALAR 
MONOPOLE COLLECTIVE MOTION 



Collective motions are easily identified in time- 
dependent approaches as periodic oscillations of the as- 
sociated collective variables. For L = excitations, an 
initial perturbation preserving radial symmetry induces 
a collective vibration of the nucleus with specific frequen- 
cies that can be extracted from a Fourier analysis of the 
oscillation pattern. The collective vibration of a spherical 
nucleus in small amplitude is typically associated with a 
unique frequency, the giant monopole resonance, whose 
energy is directly related to the second derivative of the 
energy functional with respect to the root-mean-squarc 
radius, that is, to the nucleus incompressibility coefficient 
1- 

The advantage of explicit time-dependent approaches 
compared to the more standard RPA analysis is that 
they allow multiple excitations of giant resonance and 
coupling to more complex high-energy channels [HI, [I?} ■ 
Since at large amplitude the monopole oscillation evolves 
toward radial expansion followed by nuclear multifrag- 
mentation, we may expect a correlation between the 
oscillation frequencies and the expansion speed which 
may affect the fragmentation pattern (6lj . In the stan- 
dard time-dependent treatment of collective excitations 
IH H H3, the initial state = 0)) is prepared by 
applying the generator of the collective excitation un- 
der study, Uo(k) = cxp(ifc<5 n ), on the ground-state wave 
function (fx, . . . , r^O) , where Qo = ? 2 for the isoscalar 
monopole. In this case the Fourier transform of the 
temporal variation of the collective response A a/ (f 2 ) = 



FT[r](w) 



dt AV \f 2 ){t)e 



(9) 



tional to the strength function 

3 (FT[r](w)) ex S(w) = ]T |<0|f»| 2 <5(^ - E n ) , 

71 

(10) 

meaning that the excitation probability of the different 
frequencies is a direct measure of the transition proba- 
bility to the excited states \n) of energy E n , induced by 
the excitation operator Qq. 

This treatment cannot be applied to the AMD model. 
Indeed the application of the monopole operator Uo(k) 
to the single-particle wave function, Eq. produces a 
new Gaussian state with centroid K and width a': 



~, _ l + 2ika* t 

1 ~ TTlFkP 



1 + 2ika* 



1 



(11) 



which, because of the nonzero imaginary part of the 
width, docs not belong to the AMD variational space. 
This simple observation already suggests that the dynam- 
ical evolution of the width may play a significant role in 
the breathing mode. 

However, we can still study compressional collective 
modes in AMD by choosing the initial state as a com- 
mon translation in the radial direction of the Gaussian 
centroids with respect to their ground-state value, 



R(&(t = 0)) = (1 + A)K$) 



3 (&;(* = 0)1 =9f$). 



(12) 



in the small- amplitude limit (k — > 0) is directly propor- 



This initial state, as any other modifying the root- 
mean-square radius of the nucleus, will excite the same 
monopole eigenmode as the application of the monopole 
boost Uo(k) |33j ] . This is however not true as far as the 
spectral function FT[r](w) is concerned. Indeed in the 
case of the AMD model where the Gaussian centroids 
are the only variational parameters, this initial state cor- 
responds to minimizing the total energy in the presence 
of a perturbative external field — Xr 2 for some chosen val- 
ues of A. It is possible to show [H, [62| that including an 
external harmonic field V ex t = \r 2 9(—t) in the TDHF 
dynamics is indeed an alternative way to prepare the ini- 
tial monopole state. Such an initial condition produces 
in the linear response regime a spectral function differ- 
ent from Eq. (|10p . where the probability of the different 
collective states is weighted by the inverse of the state 
energy. 

For this reason we will concentrate on the energy of 
the collective states, which appears to be independent 
of the chosen initial condition, rather than on the de- 
tailed shape of the strength function. Another argument 
for that is given by the fact that a realistic description of 
the strength function and the associated energy weighted 
sum rule (EWSR) requires a proper treatment of the 
damping of collective motion, which in turn demands the 
inclusion of higher order correlations in the form of col- 
lisions and branching [48| . Finally, since the equations 
of motion of AMD and FMD are obtained by restricting 
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the TDHF single-particle wave functions to Gaussians, 
the strength function obtained by AMD or FMD will 
approximately satisfy the corresponding RPA sum rule 
only if the restriction on the single-particle wave func- 
tions is not critical. All these difficulties are much less 
critical as far as the frequencies are concerned. For all 
the cases presented in this paper, the time evolution of 
the root-mean-square nuclear radius presents only a few 
well-distinct frequencies, meaning that the collective en- 
ergies can always be unambiguously recognized. 




20 



5 10 15 20 

E* [MeV] 

FIG. 1: (Color online) AMD prediction for the frequency fkd 
of the monopole oscillation as a function of the oscillation 
energy for 12 C, 16 0, and 40 Ca with the SLy4, SIII, and Gogny 
interactions. 

In AMD calculations for monopole oscillations of 12 C, 
16 0, and 40 Ca, the spectral function FT[r](w) shows only 
a single peak and thus a single frequency in the oscilla- 
tion. The peak frequency huj as a function of oscilla- 
tion excitation energy, which is directly related to the 
oscillation amplitude, with the three different interac- 
tions (SLy4, SIII and Gogny) are shown in Fig. [T] A 
global independence of the peak frequency on the exci- 
tation energy is observed in the case of 40 Ca, whereas 
the peak frequency gradually decreases for 12 C and 16 
at high excitation energies (^5 MeV). For low excitation 
energy (small-amplitude) oscillations (j5 MeV), the peak 
frequency varies abruptly as a function of excitation en- 
ergy for the 12 C case if the Gogny interaction is used. A 
clear, even if much less spectacular, frequency change for 
small amplitude oscillations is also observed for the other 
effective interactions and the other nuclei. This behavior 
will be discussed in detail in the next section. However, 
this behavior at small amplitudes is not likely to strongly 
affect the energy of one-phonon states when the collec- 
tive motion is quantized. The physical frequency that is 
relevant to the breathing mode is rather the one at the 



excitation energy of the mode itself, since the percentage 
of the one-phonon state is maximized at this energy. 

For the case of 40 Ca, different theoretical studies ex- 
ist on the GMR in this nucleus [H, [64[ , which can thus 
be considered as a benchmark for our calculations. The 
ground state of 40 Ca is obtained in AMD as a structure- 
less spherical distribution of neutron and proton densi- 
ties. The resulting monopole peak is found at an en- 
ergy hu> in reasonably good ag reement with TDHF (32j 
(hu: = 22.1 MeV) and RPA [§4( (ftjj = 21.5 MeV) calcu- 
lations using the same SLy4 effective interaction. 

The SIII interaction is employed here for the purpose 
of exploring the sensitivity of the monopole response to 
the nuclear matter properties. As expected, as the in- 
comprcssibility A'oo increases from 230 MeV (SLy4) to 
355 MeV (SIII), all the peak frequencies increase, show- 
ing that indeed collective compression and expansion can 
be addressed in molecular-dynamics-based formalisms. 

The peaks obtained by AMD turn out to be without a 
physical width (5- function-like) . In principle the energy 
of the modes is above the threshold of particle emission 
and the motion should be damped by the coupling to 
the continuum [H, [28|, [32[ . This is however not observed 
in the simulation, showing that the AMD wave function 
ansatz is not flexible enough to properly describe particle 
evaporation. The inclusion of higher order correlations 
is also required for a proper treatment of the damping of 
collective motions. In future studies it would be interest- 
ing to investigate whether the AMD models are able to 
reproduce the width of the collective resonances by intro- 
ducing the two-body collisions and quantum branching 



V. DIFFERENT MODES AND THE ROLE OF 
THE WIDTH 

In AMD calculations, the inverse width parameter v 
is an external parameter and the calculations shown in 
Sec. II VI were performed using the optimal value f g . s . (ta- 
ble [Hj , namely, the one that gives a minimum ground- 
state energy. However, there is no strong reason to adopt 
i/g. s . for the oscillation calculation if the collective ex- 
cited states have considerably different structure from the 
ground state. For example, if the collective excited states 
have a strong clustered structure, the optimal choice of 
v for individual clusters may be more suitable than v^ a , . 
Furthermore, in some applications of AMD such as to 
nuclear matter and heavy-ion reactions, it is not practi- 
cal to adjust v for individual nuclei that appear in the 
system. For these reasons, it is important to understand 
how the characters of monopole oscillation depend or do 
not depend on the choice of which we will study in this 
section. We take the 12 C calculation as a reference in the 
following, but qualitatively similar results are observed 
for 16 and 24 Mg. 

Figure [5] shows the position of the peak frequency in 
the spectral function FT[r](w) of 12 C by AMD with SLy4 
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FIG. 2: (Color online) AMD prediction for the peak posi- 
tion hui of the monopole oscillation for 12 C with the SLy4 
interaction as a function of the oscillation amplitude for dif- 
ferent values of the inverse width parameter v (fm -2 ). The 
excitation energies corresponding to oscillation amplitudes in 
the case of the optimal width v = 0.155 fm -2 (full line) are 
indicated by arrows. 



as a function of oscillation amplitude for different inverse 
width parameters v. The full line shows the result of the 
ground state optimal width {y^. s . = 0.155 fm -2 ). The ex- 
citation energies corresponding to oscillation amplitudes 
for this line are also indicated in the figure. 

In the region of the amplitude greater than 0.4 fm, 
different choices of v give almost identical oscillation fre- 
quencies, and in particular the prediction of GMR en- 
ergy is almost independent of the choice of v. The fre- 
quency gradually decreases as the amplitude increases. 
This behavior can be physically understood, since above 
the breakup threshold the expansion can be followed by 
a breakup of the system and not by a recompression, cor- 
responding to an oscillation period that tends to diverge 
(that is, a frequency that tends to zero). 

In small-amplitude oscillations (<0.4 fm), the oscilla- 
tion frequency varies abruptly as a function of the am- 
plitude for some choices of v. A very peculiar behav- 
ior is observed in particular for the case v = 0.17 fm~ 2 
in the small-amplitude regime. With this choice of the 
width parameter the monopole frequency shows a well- 
pronounced minimum at a specific value of the oscillation 
amplitude (0.04 fm). In the following we will show that 
this is due to the presence of two possible distinct oscilla- 
tion regimes in the AMD dynamics by investigating the 
motion of wave packet centroids. 

The motion of Gaussian centroids gives some useful 
information on the time evolution of the system, though 
the centroids do not correspond to the physical positions 
of nucleons owing to antisymmctrization. In the AMD 
dynamics of 12 C, three ensembles of four Gaussians with 
different spin-isospin (p f, p n *f, and n \.) take al- 
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FIG. 3: (Color online) AMD density distributions of 12 C with 
the SLy4 interaction projected on the plane perpendicular to 
the symmetry axis (x-y plane) for (a) the ground state and 
the maximum amplitude stage of the monopole vibration at 
(b) E* =0.1 MeV and at (c) E* = 5 MeV. The cross marks 
in the figures indicate the parameter positions of the effective 
"a"s (See text). 



most the same time-dependent value for their b parame- 
ter. The global 12 C configuration shows therefore a "3-a- 
triangular" structure during the time evolution. The den- 
sity distribution does not show spatial clustered structure 
at very small amplitude [Fig. E)Jb)] because of the large 
overlap among the "a"s and the antisymmetrization ef- 
fect among the nucleons, but these structures are clearly 
visible for larger amplitude vibrations [Fig. E[c)]. We 
will speak of effective "a"s in the following, where the 
word "effective" means that the associated four nucleons 
do not form three independent a particles because of the 
global antisymmetrization of the wave function. 




v=0.17fm" 3 



R a (fm) 



FIG. 4: Phase-space trajectory of the AMD monopole oscil- 
lation of 12 C from the evolution of the variational parameters 
of one "a" (see text) for (a) v = 0.13 fm" 2 , (b) 0.16 fm" 2 , 
and (c) 0.17 fm -2 with SLy4. Each contour corresponds to a 
different initial expansion, the largest circle being associated 
with the largest amplitude oscillation. The global oscillation 



pattern of 
(see text). 



'C can be deduced from symmetry considerations 



The active degrees of freedom in this motion are the 
variational parameters of one of these effective "a"s in 
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C, namely, the Gaussian centroid parameters 

1 1 4 

i? a = i ($ a |r1d> Q ) = I ^3?(&, 

i=i 

4 



(13) 



the other two effective "a"s move in phase to keep the 
triangular symmetry owing to the symmetry of the exci- 
tation. R a — R a u r and P a = P a u r oscillate along the ra- 
dial direction u r = R a (t = 0)/\R a (t = 0)|. The trajecto- 
ries of R a and P a for three different choices of the inverse 
width are shown in Fig. 2J The ground-state (R a ,P a ) 
values are indicated by the crosses in Fig. [T] In the case 
of a wide Gaussian width ^ = 0.13 fuG 2 [Fig. 0|a)] the 
position of the effective "a" s in the ground state coincides 
with the center of mass of the nucleus. The trajectories 
circulate in phase space around (R,P) — (0,0) and the 
other two effective "a"s move in phase so that this mo- 
tion corresponds to the crossing of "3-a" as schematically 
shown by the cartoon of Fig. [SJa) . For narrower Gaus- 



as in Fig. |4|c), the excitation from the ground state, 
which has a finite "a" clustering, should couple to the 
excitation from the parity-reflected ground state. The 
lack of such a coupling in the present calculation may 
be a reason of the low oscillation frequency around the 
transitional situation between Figs.JSJa) andJSJb). How- 
ever, if the oscillation amplitude is large (namely, if the 
"deformation" of the intrinsic state I^amd [<?(*)]) is large 
enough), we can expect that the effect of the parity and 
angular momentum projection is not so important. It 
should also be noted that a kind of many-particle many- 
hole excitations (namely, cluster excitation) is included 
in the one-phonon state of the oscillation if it is extracted 
from the time-dependent many-body state I^amd [?(*)]) 
or P(0 + )|*amd [?(*)]) around the GMR amplitude. 
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35 
30 
>25 

CD 

1,20 
|15 
10 
5 




(a) 


(b) 12 CSLy4 







FT [arb. units] 



5 10 15 

E [MeV] 



FIG. 5: Schematic representation of the two possible oscilla- 
tion regimes observed in AMD calculations varying the Gaus- 
sian width and/or the motion amplitude. 

sian width v > 0.16 fnG 2 , the ground state R a becomes 
finite, indicating an explicit "a" structure even in the 
ground state. In this case the small-amplitude oscilla- 
tions [Fig. HKc)] correspond to vibrations around this fi- 
nite value, leading to a trajectory confined to the positive 
i? Q -coordinate space. This motion can be visualized as 
three effective "a"s oscillating around their ground state 
values as schematically indicated by Fig. EJb). Increas- 
ing the excitation energy, and therefore increasing the 
amplitude of the motion, we find trajectories circulating 
around (R a , P a ) = (0, 0) as in the case of v = 0.13 fnG 2 . 
The frequency minimum for v = 0.17 fnG 2 in Fig. [2] is 
observed at the transition between these different oscilla- 
tion regimes. The small-amplitude oscillations observed 
with v = 0.16 fnG 2 [Fig. H£b)] are located at the transi- 
tion between the two motions, which explains the strong 
v dependence of the peak frequency for the small ampli- 
tude oscillations in Fig. [2] 

As the time-dependent state |^amd[<?(*)]) violates 
the symmetries, a more desirable approach is to ap- 
ply the time-dependent variational principle to the state 
P(0 + )| x f , AMD['?]) after the parity and angular momen- 
tum projection. For the oscillations with two regimes 



FIG. 6: (a) The spectral function FT[r](u) in arbitrary units 
of the FMD monopole oscillation in small amplitude for 12 C 
with the SLy4 interaction, obtained by the same initial state 
as in the AMD calculation, (b) The peak position Tiw as a 
function of the oscillation amplitude. 

From these discussions, it is very interesting to see the 
results of FMD, which explores dynamically the different 
values of wave packet widths {a} as well as the centroids 
{b}. Figure [6] shows the position of the peak frequency 
hw in the spectral function FT[r](w) of 12 C by FMD as 
a function of oscillation energy, where the initial state 
is prepared by Eq. (fT2|) . The spectral function in small 
amplitude is shown on the left side of this figure. The 
widths of the peaks are only a numerical artifact coming 
from the time filter (28|, which is applied to compensate 
for the finite integration time in the Fourier transform. 
In the case of FMD, there are two peaks in the spectral 
function of 12 C and the two lines in Fig. |5|b) correspond 
to the frequencies of these peaks. The behavior of the 
lower frequency follows closely the behavior of the one 
observed in AMD in the excitation energy regime where 
no sensitivity to the width was observed. However, a 
smooth behavior is also obtained at the lowest ampli- 
tude, showing that indeed two types of vibrations could 
be physically present. 

For comparison, the experimental monopole strength 
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distribution Q shows a fragmented structure composed 
of at least four different peaks between 13 and 45 MeV of 
excitation energy exhausting around 27% of the EWSR. 
The average energy in this excitation energy range is 
mi /mo = 21.9 MeV with an root-mean-square width of 
4.8 MeV, in qualitative but not quantitative agreement 
with our results. The experimental data also show a low- 
energy peak that is associated with the excitation of the 
Hoyle state. No strength at these low energies is observed 
in our calculation, consistent with the fact that effects 
beyond mean field have been shown to be necessary to 
correctly reproduce the Hoyle state [H, 5s| • 

To understand the physical meaning of the two peaks 
observed in the calculation, let us investigate the dynam- 
ics of variational parameters of 12 C in the model. Similar 
to the case of AMD, three ensembles of four Gaussians 
with different spin-isospin (p ^, p 4-, n t, and n \.) take 
almost the same time-dependent values for their b param- 
eters and constitute a "3-a-triangular" structure. The a 
parameters associated with each of the three "a"s oscil- 
late in the same way. Therefore, the active degrees of 
freedom in the FMD motion are given by the centroid of 
the "a"s like in AMD, 



It 



9(04)9(6. 



3? (Oj) 



4 3(6, 

i=i v i > 

together with the width of the effective "a"s, 



(14) 



°R = ($ a \r 2 \$c 



E 



3|a, 



o P = ^ a \f\^ a ) = J2 



^25i (t»i) 
3fi 2 



(15) 



Figure [7] shows the oscillation frequencies associated 
with these different observables. We can see that the 
first peak present in the monopole response (Fig. [5]) is 
associated with the time variation of the centroid of "a" 
R a . This physically corresponds to an oscillation of the 
positions of three "a" clusters around the center of the 
nucleus [see Fig. @Ja) and Fig. [5ja)]. In the absence of 
the width degree of freedom, this is the only collective 
monopole motion accessible to the system, and this is 
why the AMD response contains a single peak. The pre- 
diction of AMD and FMD concerning this peak is in good 
agreement because it comes out that 3(a) ~ during the 
whole dynamical evolution. 

The second peak at higher energy observed in the 
monopole response (Fig. [5]) is associated with the time 
variation of the widths. It is physically associated with 
the coherent vibration of the constituting "a" clusters, 
as shown by the fact that this frequency appears in the 
Fourier transform <jR(t). This mode is very similar to the 




15 20 25 
fico [MeV] 



40 



FIG. 7: (Color online) Fourier transforms in arbitrary units of 
FMD monopole oscillations ( 12 C with the SLy4 interaction) 
associated with the square of the position of the effective "a" 
from the 12 C center of mass i? 2 (t) (dashed line); the coor- 
dinate width of the effective "a" cm(t) (thick line); the dif- 
ferential time evolution of the proton and neutron for square 
radius (i? pro ton - -Rncutron) (thin line). 



standard breathing mode of medium-heavy and heavy 
nuclei, where all nucleons vibrate coherently. 

In conclusion, there arc two possible isoscalar 
monopole oscillations in the FMD dynamics of an excited 
12 C nucleus: (1) cluster moving and (2) cluster breathing. 
This result is not unique for 12 C, but it is also found for 
16 O. The cluster-originated vibration modes are system- 
atically found for all light nuclei that in the framework of 
the AMD and FMD model exhibit a pronounced cluster 
component. Vibration (1) finds its natural prolongation 
at high amplitude in the phenomenon of direct cluster 
breakup, which is systematically observed in our calcu- 
lation at very high excitation energy out of the linear 
response regime. 



hw (MeV) 


initial state (1111) 


initial state (1121) 


SLy4 


SIII 


SLy4 


SIII 


12 C 


24.2 


30.4 


20.2 & 24.2 


20.8 & 30.7 


16 


26.1 


32.4 


26.2 & 26.8 


26.8 & 32.8 


40 Ca 


22.1 


27.6 


20.9 


27.5 



TABLE III: Energy of the collective resonances fiu recog- 
nized as peaks of the spectral function FT[r](oj) obtained by 
FMD for different selected light nuclei ( 12 C, 16 0, and 40 Ca) 
and effective interactions (SLy4, SIII, and Gogny) in the low- 
amplitude limit with two different initial states. 

Table IIIII shows the summary of the peak frequen- 
cies hui in the spectral function FT[r](w) obtained by 
FMD for different selected nuclei and effective interac- 
tions in the small-amplitude limit with different initial 
excitations. We can see that both the cluster moving 
and cluster breathing modes are excited when the ini- 
tial state is prepared by a radial translation of the wave 
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packets' centroids [Eq. (|T2|) ]. whereas only the cluster 
breathing mode is observed if the standard monopole 
boost [Eq. (fTTj) ] is applied to the ground state FMD wave 
function in the small-amplitude limit. This suggests that 
the cluster moving mode in 12 C may be due to a cou- 
pling between the initial excitation [Eq. (|12j) ] applied on 
a deformed ground state and the monopole mode built 
on top of it, owing to the nonlincarity of the FMD dy- 
namics at finite amplitude. A similar effect was already 
reported in TDHF calculations in Ref. (33|. Moreover, 
in the case of 12 C, the experimental threshold energy for 
3-a breakup is given by E* = 3E g . s .(a)- E g . s .( 12 C) = 7.3 
MeV which is below the energy associated with the col- 
lective frequency of the cluster moving mode of AMD 
and FMD, even considering the Coulomb barrier, which 
is expected to be around 2 MeV [l3| . The observed vibra- 
tion appears then more as a virtual diffusion state than 
as a resonant state. This is however not the case for all 
cluster moving modes. As an example we show in Fig. [5] 
the Fourier transform FT[r](w) of the FMD oscillation in 
small amplitude for 24 Mg. 



24 Mg SLy4 
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tative discussion. 

For comparison, the experimental monopole strength 
distribution [ll[ of this nucleus shows a very broad struc- 
ture peaked at around 21 MeV with a root-mean-square 
width of around 5-7 MeV, exhausting around 80% of the 
EWSR. Here again, only qualitative conclusions can be 
advanced. If the complicated structure of the spectral 
functions that are systematically observed in light nu- 
clei are certainly an indication of a global loss 
of collectivity, our study suggests that a multiple-peak 
monopole spectrum can also be an indication of clus- 
tered states. A prolongation of this work requires the 
inclusion of collisions and branching in the molecular dy- 
namics [48| , which would allow a quantitative prediction 
of the strength function as a possible observable to evi- 
dence cluster structures. 

Other collective modes can also be investigated with 
our approach by choosing an appropriate initial state and 
analyzing the associated observable. In fact, in the oscil- 
lation calculation of 12 C discussed here, when adopting 
the initial state (p~2|) . we also find another peak around 
33 MeV by looking at the Fourier transform of the width 
variable of a proton wave packet. This extra peak is asso- 
ciated with the isovector monopole mode, which is clearly 
seen in the Fourier transform of the differential time evo- 
lution of the proton and neutron root-mean-square ra- 
dius (-Rproton - -Rncutron) (the thin line in Fig. [7]). In- 
deed the chosen initial state, where all the wave packets 
have the same width, contains some isovector excitation 
in the FMD model because protons and neutrons have 
slightly different widths in the ground state owing to the 
Coulomb force. 



VI. SUMMARY 



FIG. 8: The spectral function FT[r](oj) in arbitrary units of 
the FMD monopole oscillation in small amplitude for 24 Mg 
with the SLy4 interaction. 

The ground state of 24 Mg is found to be a deformed 
prolate shape and it is natural to find a fragmented 
structure in the spectrum function of such a nucleus 
since the monopole oscillation is strongly coupled to the 
quadrupole oscillation. However, the analysis of the vari- 
ational parameters {a, b} reveals that once again the 
peaks are related to different motions of the underly- 
ing quasi-clustcrs. Specifically, the origin of the three 
major peaks (hui = 16, 22, and 25 MeV) are seen to 
correspond with 3 "Be", 2 "Be" + 2 "a", and 6 "a", re- 
spectively, where "Be" denotes the effective cluster for 
8 Be as "a" for 4 He. These peak frequencies are be- 
low the experimental threshold energy for 6-cv breakup, 
E* = 6E g . s .(a) - E g . s .( 24 Mg) = 28.5 MeV and therefore 
these resonances may physically exist, although the im- 
plementation of the spin-orbit interaction and the projec- 
tion over angular momentum and parity before the vari- 
ational procedure will be necessary for a further quanti- 



In this paper we have studied collective isoscalar 
monopole oscillations for some selected light nuclei in 
the framework of the AMD and FMD time-dependent 
approaches. 

The study of GMR in 40 Ca reveals that for medium- 
sized nuclei molecular dynamics models are in good 
agreement with TDHF calculations. A well-defined col- 
lective breathing mode is observed at a frequency that 
depends on the incompressibility modulus of the under- 
lying effective interaction. 

Clear collective modes are observed also for nuclei as 
light as 12 C, but the collective modes in light nuclei ap- 
pear dominated by quasimolecular cluster structures. We 
have investigated in great detail the motion of 12 C as an 
example. The 12 C isoscalar monopole dynamics leads 
to a spectral function containing two distinct frequencies 
associated, respectively, with an oscillation of the effec- 
tive three "a" constituents around the nucleus center of 
mass and with a coherent vibration of the underlying 
clusters. This latter frequency sensitively increases with 
increasing incompressibility of the associated equation of 
state, but does not follow the empirical oc^l -1 / 3 law of 
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the well-known breathing mode in medium-heavy and 
heavy nuclei. Only the first mode is accessible to AMD 
calculations, whereas the inclusion of the width degree 
of freedom appears essential to describe the second one. 
More complicated spectral functions can be observed in 
our calculations depending on the cluster structure of the 
nucleus. The presented calculations do not contain two- 
body collisions nor quantum branching, meaning that 
they cannot describe the decay of the observed collective 
motions, which are expected to considerably increase the 
fragmentation of the monopole strength. However, the 
dominant modes should be preserved by such improved 
calculations [65[. 

Molecular dynamics calculation such as AMD and 
FMD therefore will be a powerful tool to study the inter- 



play among clustering structures, collective resonances, 
and multiple breakup of various nuclei in the nuclear 
chart including unstable nuclei. 
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